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1. Introduction 

This paper continues the study of static self-gravitating elastic bodies in Einstein 
gravity which was initiated in fP. In that paper, connected bodies were considered. 
Here, we shall consider the problem of constructing static, elastic, many-body sys- 
tems in Einstein gravity. 

In Newtonian gravity, one may easily construct static many-body systems con- 
sisting of rigid bodies. Newton showed that the potential inside a spherical, ho- 
mogenous shell of matter is constant. It follows that a body placed at rest anywhere 
inside the shell, will remain at rest. This example generalizes to hollow ellipsoids, 

The problem of constructing static, self-gravitating many-body configurations 
consisting of elastic bodies in Newtonian gravity was studied in [3]. One of the 
results proved there is that under certain conditions it is possible to place a small 
elastic body at a stationary point of the Newtonian potential of another body. We 
remark here that the above mentioned example of a body at rest inside a spherical 
shell is only possible in the case of elastic bodies if the smaller body is spherically 
symmetric and placed at the center of the spherical shell, see section l4?3l below. 

The two main steps in the construction used in [1] are the application of the 
implicit function theorem to construct solutions of a certain projected version of 
the reduced, coupled system of Einstein and elastic equations, and an equilibration 
argument which showed that the solutions thus constructed are solutions of the full 
system. The implicit function theorem is used to deform from a relaxed body with- 
out gravity, i.e. with Newton's constant G = 0, to a solution of the above mentioned 
projected, reduced system with G ^ 0. As was observed in p] the linearization of 
the elastic system at a relaxed configuration has a cokernel corresponding precisely 
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to the Killing fields of a Euclidean metric. It is in order to cancel this degeneracy 
that one is led to applying a projection to the system. 

In generalizing the method used in the case of a connected body to the case of 
multiple bodies, it turns out that the linearized elastic system now has a cokernel 
corresponding to the product of the cokernels corresponding to each body. However, 
the equilibration argument developed in [T], which is essential to show that the 
solution constructed is a solution of the full coupled system of Einstein and elasticity 
equations, requires that one has a solution corresponding to a system which is 
projected along only one of the factors of this cokernel. 

In order to achieve this situation one must in the case of multiple bodies, after 
solving the system resulting from projecting out the cokernel, perform a partial 
equilibration making use of the additional degrees of freedom one has in a multi- 
body system provided by the relative position and orientation of the component 
bodies. This step requires certain non-degeneracy conditions on the Newtonian 
potential of the bodies, which are completely analogous to the ones considered in 
i- 

After the partial equilibration, the only remaining degeneracy corresponds to 
the cokernel associated with one of the component bodies, and the equilibration 
argument from [T] can be applied to finish the argument. 

In sections we present the field equations for self-gravitating elastic bodies 
and prove the existence of solutions to this system, for disconnected bodies, given 
certain non-degeneracy conditions on the relaxed reference system. In section |4l 
we give a few concrete examples where this construction applies. 

We remark that for simplicity, the results in sections [2H1] are stated in terms 
of self-gravitating two-body systems. However, the method and results apply with 
essentially no changes to the construction of n-body systems also for n > 2. The 
point which needs to be noted here is that instead of partially equilibrating one of 
the bodies, cf. section [3l for the case of a system of n bodies, it is necessary to 
partially equilibrate n — 1 of the bodies, by a procedure completely analogous to 
the one presented here. We leave the details to the reader. 

The construction of static solutions of the Einstein equations representing discon- 
nected elastic bodies raises the question of the necessary conditions for the existence 
of such configurations. In particular, what condition rules out the existence of two 
static bodies in Einstein gravity? 

In Newton's theory of gravity it may seem naively clear that two bodies can not 
remain at rest under their mutual gravitational force. However, it is important to 
note that in order to be able to make such a statement, one needs a separation 
condition. For example, consider a hollow sphere and place a small spherical body 
in the center. As is well known, the Newton potential is constant in the interior of 
a sphere, and therefore such a configuration is static. 

A natural condition which rules out the existence of a static two-body system 
Bi U B2 in Newtonian gravity is given by the existence of a plane which separates 
the two bodies Bi and B2- In this case it is easy to see that no static solution can 
exist. The self force of Bi vanishes. Consider now the component of the total force 
generated by B2 acting on Bi, which is orthogonal to the separating plane. This 
component is clear non-vanishing and hence the bodies must begin to move. 

Now consider the same question in Einstein's theory of gravity. As we show in 
this paper, there are static two-body configurations. However, in Einstein gravity. 
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it is not clear how to generalize the separation condition discussed above from the 
Newtonian case. In the paper [4] , it was proved that there are no static solutions of 
Einstein gravity which admit a separating, complete, totally geodesic hypersurface 
E. Here, separating can be taken to mean that S lies in the complement of the 
bodies. In particular, this rules out a static two-body configuration where the 
bodies are separated by a totally geodesic hypersurface. 

The problem under discussion, which is rather simple in Newton's theory, be- 
comes a challenge in the context of Einstein's theory of gravity. One explanation 
for this additional difhculty is that in Einstein gravity, we loose the concept of force. 
Furthermore, the self field is no longer a useful concept because the field equations 
are nonlinear. Miiller zum Hagen wrote a PhD thesis on this problem and has a 
result for two separated axially symmetric fluid bodies, cf. 

2. General setup and solution of the projected system 

We adopt the notations and conventions of \Xs. We here consider the situation 
where the reference body has two connected components. 

The bodies are assumed to be disjoint domains in Rg, the extended body, and 
to have free boundaries. In particular, we consider, as in [1], the zero traction 
boundary conditions. We consider configurations / : ]R| ^ and deformations 
(p : R| as in the case of a connected body. Here Rjl is the space manifold. 

The same constitutive relations as in [Tj are assumed to hold. 

We will assume that B satisfies the Newtonian equilibrium condition 

f £_'^^v = 0, £=1,2 (2.1a) 

Av = inlxB (2-lb) 

for any Killing field ^ on R^, see the discussion in section l3^ for explanation of this 
condition. Here e is the rest mass term in the relativistic stored energy function, 
see [11 section 3.3] for discussion. 

The field equations are the same as in the case of a connected body, namely in 



Lagrangian frame 

V,ie^a^) = e^{ne-ai')V.,U in f-\B), a,^ n,\j-i^oB) ^ (2.2a) 

AhU = A7:Ge^ine~a/)xf-iiB) in R| (2.2b) 

= 8^TG{e^J - e^ay Xf-HB)) in R| (2.2c) 

and in material frame 

V^(e^a/) = e^[e-^]9^ mB, s/n^laf? = (2.3a) 
n 

Ajj = AirGe^ {ne - ai^)xi3 in Rg (2.3b) 

-^Aj^+Q,,(h,dh) = 2(V,C/)(VjC/) - 87rGe^((Ty ~h~ai')xi3- (2.3c) 
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2.1. Analytical setting. Let Bi = W^-p{B) x M^f ^ x Ef^, and let B2 = [LP{B) x 
B^^^^P'P{dB)] X £^_2 ^ '^?-2- Then, Si is a Banach manifold, and B2 is a Banach 
space. 

The residuals of equations (|2.3p define a map : R x Si ^ ,62, ^ = -^(G, Z), 
where we use Z = {(j),U,hij) to denote a general element of Bi. We assume that 
is a diffcomorphism onto its image. Thus, has components = {J-^,J^u,J^h), 
corresponding to the components of B2, given by 

^0= (^VA(e%^)-e^[e--^]9^, tTaeia^-^hA^ (2.4a) 

^^Jl-47rGe^{fie-a/)XB (2.4b) 
= - ^ AJ^ + Qij (h, M) - 2V~UV~U + SttGc^ (a.y - 7^ a; ' )xe (2.4c) 

Recall, cf. [l] section 3.2], that in setting up the field equations in the material 
frame, an extension of from B to is used, cf. 1, Lemma 4.1]. The argument 
used to construct (j) carries over without change to the case of a body with several 
connected components. 

The equation to be solved is J-{G, Z) = 0. The material form of the reference 
state is given by 

Zo = (i,0,(5y oi) e Si. 

The map T defined by ()2.4p is easily verified to satisfy J^(0, Zq) = and to map 
Bi B2 locally near the reference state Zq. As discussed in [T], the map T is 
smooth. 

The Frechet derivative 02^(0, Zq) can be represented as the matrix of operators 

DrpJ^rj, DuT^ DhJ-^A 
A 

-ia; 

(where the entries are evaluated at {0,Zq)). In particular, the matrix is up- 
per triangular, and the diagonal entries are isomorphisms, with the exception for 
D^J-^{Q, Zq) which is Fredholm with nontrivial kernel and cokernel. 

Let (jje denote the restriction of 4> to Bi. Then we have in a natural way J- = 
{T^^,T^.2,J^u,^h), and £'0^^0(0,^0) takes the form 

fD^.T^, \ 

Let 5a denote any combination of the Frechet derivatives of (Tj^, evaluated at 
(0,Zo). Assuming we use a coordinate system X"^ where V123 — 1, we have the 
relations 

0=/ CdAiSa,^)- f ei6a,^)nA, e ^ 1,2 

JBt JdBe 

where Bi are the connected components of B and is the outward normal. This 
can be interpreted as saying that due to the natural boundary conditions, the 
linearized elasticity operator, restricted to each component of B is automatically 
equilibrated at the reference configuration {0,Zo). It follows that the cokernel of 
the operator 

D^T40,Zo) : W^ P{B) [LP{B) x B^-^^P^P{dB)] 
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consists of the space 

3^1® 3^2 

where for I — 1,2, yi is the space of Killing fields on Bi considered as a subset of 
(Mg, (^e), where = i*i5 is the Euclidean metric on Mg induced from the Euclidean 
metric b on R| . The fact that the range and cokernel of JF^ consists of fields on B 
is due to the fact that we defined T by passing to the material frame. 
Similarly, the kernel of D^J-'^(0, Zq) is the space 

where is the space of Killing fields on Bg, considered as a subset of (Rg, Sis)- 
Hence, in view of ellipticity, the operator 

D^T^{0,Zo) : W^'P{B) -> [LP{B) x B^-^/P^P{dB)] 

is Fredholm with the finite dimensional kernel and cokernel discussed above. 

2.2. Projections. Introduce the projection operators Pg^ : B2 B2, which acts 
as the identity in the second and third components of B2 and is defined in the 
first component of B2 as the unique projection along the space of Killing fields on 
{Bi,Sis), onto the range of D^^J-'^^{0, Zq), which leaves the boundary data in the 
first component of B2 unchanged. If we consider each component Bi separetely, 
the situation is analogous to the one for the case of a connected body, so that 
lPi3f£'0<.-^0f (0, Zq) is a surjection. 

We now go back to letting the projection operator act on all components of 
and define V131I1B2 = Pfii © ^62 by combining the projections associated with each 
component of the body. Then we have that IPbiuB2 ^2-^(0, Zq) is a surjection. 

The projected system has the property, as in the case of a connected body, cf. 
[1] that, 

PBiU62^2^(0,^o) 

is a surjection with finite dimensional kernel. For this reason, the implicit function 
theorem can be applied more or less directly to construct solutions to the first 
system. 

The data for an Euclidean motion is given hy A = (a*, (3ij), where a' e M'^ is a 
translation vector and /J^ is an orthogonal matrix. The motion A acts in M.g by 
X 1-^ P{x + a). Denote the group of Euclidean motions by A, and let / e ^ be the 
identity. We shall consider situations where to leading order cf) maps B2 to v4i(Z?2)- 
We implement this by putting conditions on the the 1-jet of </) at a point X2 G B2- 

2.3. Solving the projected equation. The following result is analogous to [U 
Proposition 4.3]. However, here we construct a family of solutions to the first 
projected system with parameters (G, A), where A A. The proof is an application 
of the implicit function theorem. 

Proposition 2.1. Let T : Bi ^ B2 be map defined by Ii2.4\ ) and let PBiue2 
defined as in section \2.'A Let Xi G Bi, X2 € B2 be given points, and let A be an 
Euclidean motion. Then, for sufficiently small values of Newton's constant G, and 
for A sufficiently close to L, there is a solution Z = Z{G, A), where Z = (</>, U , hij), 
to the reduced, projected equation for self-gravitating elastostatics given by 



PsiuB2-^(G,Z) -0, 



(2.5) 
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satisfying the conditions 

(</)-ir(Xi) = 0, S'',SciAdB]{<p-inXi) = (2.6a) 

{^-Aoiy{X2)^a\ S'',SciAdB]i<p-Aoiy{X2)=f3AB (2.6b) 

In particular, for any e > 0, there is a G > 0, such that Z — Z{G,A) satisfies the 
inequality 

\\(t> - illiy^^pCBi) + 110 - ^oi|lw2.p(g2) + ll^y - + < e- (2-7) 

3. Equilibration 

In this section we will make use of the solution to the projected system ()2.5p 
to construct solutions of the full system of Einstein equations for two static elastic 
bodies. Given a solution to the projected system (|2.5p as in Proposition 12.11 our 
first goal is to construct a family of solutions to (|2.5p which arc equilibrated on 
one component of the body. Once this is done, we are in a situation where we are 
able to apply the equilibration argument presented in [TJ section 5] to construct 
solutions to the full system of equations for the self-gravitating elastic body with 
two components. 

3.1. Partial Equilibration. For definiteness we shall focus on B2 and construct 
a curve Z = Z{G) such that 

/ Ch = o (3.1) 

for all Killing fields on M|, where 

h^VA{e^aj^)-e^[e-^]W 

n 

is the first component of J-^p. We do this by finding A ~ j4(G), with A[Q) — /, such 
that Z{G,A{G)) solves ((XT|) as an equation for A = A{G). 

3.2. The normalized force. Let (C(q,))q=i be a basis for the space of Killing fields 
on Mg. It will be convenient to solve (j3.ip by transforming to the Eulerian frame, 
using the change of variables formula as in [T], section 5.1]. 

We define the force map N = (N(„)(G, Z(G, A)))6^i, N : R x ^ M^, by 

N(„)(G, A) = / C„) ° [V,(e^a,J) - e^(ne - ai')^,U] d^ih (3.2) 
"'0(82) 

where the right hand side is evaluated at Z{G,A). The form of the force map in 
the material frame is easily found by analogy with (|3.1|) . We write this as 

We will freely make use of the form of N which is most convenient. Since there is 
a factor G in equation (j2.2bl we see that N(0, A) = 0. Hence it is convenient to 
introduce a normalied force map N for G 7^ by setting 

N = G"iN (3.3) 

We define N(0, A) by taking the limit as G 0, which is easily shown to be well 
defined, see below. It is natural to view N and N as taking values in the dual of 



STATIC MANY-BODY SYSTEMS 



7 



Lie(^) = M^. If we denote by (•, •) the pairing between Lie(^) and its dual, then 
we can write eg. 

Introducing a new potential V by GV — U , we have 

N(„) = / eL) ° [V,(e'^^G-V,^) - e«^(ne - ai')V,V] dpi^ (3.4a) 

where from (j2.2bp . V and h solve 

AhV ^ Ane^^ {ne - ai')x^(^B) in M|, (3.4b) 
G„- = 87rG(G2e[T/],, - e«^a,j X/-i(B)) in M| (3.4c) 

3.3. Newtonian Equilibrium condition. In order to evaluate N at (0, J) we 
must consider the limit limc^o N. Calculating this limit is equivalent to calculating 
the derivative 9gN(0, /). Differentiating the system (|3.4p with respect to G, at 
G = 0, we find 



N(„)(0, A) = / C [dj{Sa),^) - ed,V] , (3.5a) 

JAi(B2) 

AV = 47re(xi(Bi) + XAi{B2)) (3-5b) 

where Sai^ is a collection of derivatives of (7i^ with respect to (j). The first term in 
the right hand side of (|3.5ap vanishes identically, cf. the discussion in section [531 
see also [IJ §4.2]. Therefore we have 



N(„)(0,^)= / eic.)^^V, (3.6a) 

JAi{B2) 

AV = 47re(xi(£3i) + XAi(B2)) (3-6b) 

for Killing fields ^(q.) . We shall look for a family of solutions of N(G, A) = of the 
form A = A(G), with A{0) = /. For this to be possible, it is necessary that the 
condition N(0, /) = holds. Let Ve be the Newtonian potential of the components 
Bi. In view of the third axiom of Newton, the principle of actio est reactio, cf. [31 
§5], the self-force of a body vanishes. Applying this to each component and to the 
whole body we have 



0= / e^^V^, £=1,2, 



and 



= / f a,F2 + CVi (3.7) 



This implies that N(0, /) = takes the form 



0= / C^^Vl (3.8) 

JB2 

for Killing fields f\ We are assuming that (|3.8p holds for the reference configuration, 
cf. (|2.1ap . In particular we have 



N(„)(0,^)- / Cio.)^^Vl, (3.9a) 

JAi(B2) 

AT/i =47rexi(8i) (3-9b) 
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for Killing fields ^(a)- 

3.4. EfTect of motions on the normalized force. Next we consider the deriv- 
ative 9^N(0,/). We must consider the A derivative of equation p.9ap at G = 0. 
To do this, we must consider the effect on Ai{B2) of its motion in the potential of 
i{Bi). Let A — I + eri + 0{e^), where 77 is the infinitesimal motion with data {a, /3), 
i.e. the Killing field if ~ a* + f3'^jX^. Then we have A^^ = I — er] + O(e^), so that 

dA{A-')\^^j.v = -v 

By the change of variables formula 

/ C(a)5.^i=/ {qc.)^^V^)oA-' 

JAi{B2} Ji{B2) 

Differentiating the integral with respect to A at / in the direction rj gives 

The Lie bracket [C(q)7 ??] is again a Killing field, and hence in view of the fact that, 
by assumption, each component is equilibrated separately, cf. equation (|3.8p . the 
first term integrates to zero. This leads to 

9^N(„)(0,/).?; = - / qM„,Vn]"' 

We may view (?^N(0,/) as a linear map M^. If this is invertible, we may 

again apply the implicit function theorem and solve N = 0. 

Proposition 3.1. Assume that the reference body B C is in equilibrium in 
the sense that \2. 1 a\) holds. Let Z = Z{G,A) be the solution to PbiUB2^ — 
constructed in Proposition \2. 1\ and let N(G, A) be the normalized force map defined 
by 

Suppose that the derivative i9yiN(0, /) : — s- R® is invertible. Then there is an 
e > and a smooth map G 1— > A{G), [0, e) ^ A with 

N(G,A(G)) = 

/orGe [0,e). 

3.5. Equilibration. We are now in a position to apply the method developed in 
[1] to complete the construction of solutions to the full system (|2.3p . which then 
also gives a solution to (12. 2p . If the assumptions of Proposition 13.11 hold, then 
we may assume without loss of generality that B2 is equilibrated. Therefore we 
are in a situation which is completely analogous to that considered in [1] §5], and 
a straightforward application of the methods developed there yields the following 
result. 

Theorem 3.2. Let Z[G, A) be the solution to the reduced, projected system of 
equations for a static, elastic, self- gravitating body 

^BiiJB2^ = 0, 

constructed in Provosition \2. 1\ Assume that the normalized force map satisfies 

N(G,^(G)) = 0. 

Then in fact Z(G, ^(G)) solves the full system \2. 2\} of equations for a static, elastic, 
self- gravitating body. 
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The following is an immediate corollary of proposition 13.11 and theorem 13.21 

Corollary 3.3. Assume that the reference body B C is in equilibrium in the 
sense that i2.1a\) holds. Let Z ~ Z{G,A) be the solution to PBiue2-^ — con- 
structed in Proposition \2. 1\ and let N{G, A) be the normalized force map defined by 

Suppose that the derivative 9yiN(0, /) : ^ is invertible. Then there is an 
e > and a smooth map G ^ A{G), [0, e) A such that Z = Z{G,A{G)) is a 
solution to the full system \2/0jl of equations for a static, elastic, self- gravitating 
body. 

We have now reduced the problem of construcing a static self-gravitating two- 
body system to the question of whether the normalized force map has the property 
that 9^N(0, /) is invertible. This is clearly determined by the properties of the 
reference body B. In section S] below, we consider some particular cases. 

4. Examples 

In this section we give some examples of situations where the results developed 
in this paper apply. In view of corollary 13.31 the it suffices to that normalized 
force map of the reference body has invertible Jacobian. This condition on the 
reference body is precisely equivalent to the condition needed for the case of static 
elastic Newtonian two-body systems considered in [3^ section 5]. In general, for 
each example considered there, we have a corresponding example of a an static 
elastic self-gravitating two-body system in Einstein gravity. We shall here present 
an independent analysis of these constructions. 

4.1. Small body. Here we consider a situation analogous to the one discussed in 
[21 section 5]. Let Bi be given and let Vi be the Newtonian potential of 
Assume Vi has a non-degenerate stationary point, which we may without loss of 
generality assume to be at the origin O of the cartesian coordinate systems on 
and R|. Thus, diVi is of the form 

d,Vi = B,,x^ + Od^n (4.1) 

We may without loss of generality, after rotating the coordinate system, assume 
that Bij is diagonal, Bij — aiSij. 

We consider a test body B2. We may without loss of generality assume that B2 
has its center of mass at the origin O, i.e. 

/ i = 1,2,3 

Define 

J''^ = / x'x^ (4.2) 

We call J'^ as the tensor of inertia of B2, it should however be noted that the 
standard usage, cf. [B] section 5.3] is to define the inertia tensor as 

/ \x\^6'^ - x'x^ 
JB2 
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We will now show that there is a homothety F of M| such that Fi(S2) is equilibrated 
with respect to the Newtonian potential Vi, i.e. 



'Fi(B2 

A homothetic motion of M.^ can be written in the form 



[ e^^Vl^O, V Killing fields ^ of M|. (4.3) 

m be written in t 
Fx = XQ{x+p) 

where A > is a scale factor, Q is a rotation and p G M'^ is a translation. We 
first consider homotheties of the form F{X,p)x — X{x +p). Then by the change of 
variables formula and (|4.ip . we have for = a^, 



-Fi(B2) Ji(Bi) 

Thus, defining the normalized force map R+ x R'^ ^ R'^ by 

N,(A,p) = A-4 / d,V 

we have N(0, 0) = and 



lF(\,p)i(B2) 

5p,N,(0,0) = \B2\B,j 



Thus, if the matrix Bij is invertible, then we may apply the implicit function 
theorem to conclude that for small A > 0, there is a smooth function p{X) satisfying 
p(0) = such that with F(A,p(A))i(;B)) is equilibrated. In the following we will 
consider this case only. 

Next, we consider homotheties of the form F{X,Q) — X{Qx + p(A)). After 
applying the transformation x — > A(a; +p(A)), B2 is equilibrated with respect to 
translational Killing fields. Thus it is sufficient to consider rotational Killing fields 
5*(a;) = jX^ and motions of the form 

F(A,g) = X{Qx+p{X)). 

The change of variables formula and (14. ip gives, after taking into account the fact 
that p{X) ^ 0(A), 

f Cd;V = X^ f l3\Q^jx^B,^Q-^kx'' + 0{X'') 
JFHB2) J 

Let S0(3) be the group of rotations of R^, and consider the normalized torque map 

T : R+ X S0(3) ^ 

defined by 

T(„)(A,Q) =A-5 /■ q^^d,V,a^ 1,2,3 

J\Qi{B2) 

where C(q,)(2;) = P\a)j^'' ^ basis for the Lie algebra of S0(3), so(3) = R'^, in 
particular after raising an index we have/3^^^) =/3f3. 
We calculate T(0, Q) to be 

Ji(B2) 

where 
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expresses the fact that rotating the body by Q induces an orthogonal similarity 
transformation of the inertia tensor J. In particular, there is a Qo such that — 
QJQ^ is diagonal. Due to the fact that is skew, we have 

since by assumption B^i is diagonal. Thus, after applying a rotation to S2, we may 
without loss of generality assume that is diagonal, J*-' — p'^S"^^ so that we have 

T(0,/) = 

We now calculate 9QT(0,/)./i for fi £ so(3), i.e. /i = — /i'. We have 

^QiQJQ').^, = [^,,J] 

Let g[(3) be the space of 3 x 3 matrices, and for A G 0t(3), let Ad^ be the linear 
operator defined by AcIaB — [A,B]. Recall that the Frobenius inner product on 
the space of matrices is {A, B)jr = trAB^. We can now write 

9QT(0,/) = -|i^2|(/3,AdBAdjM>^ 

We observe that Ads maps so (3) st)m(3) for B e st)m(3). If B is diagonal, 
which is the case we are considering, then Ads maps so (3) into the three dimen- 
sional subspace of st)m(3) consisting of symmetric matrices with vanishing diagonal 
elements. Further, Ad : st)m(3) — > so(3). Thus we have Ad^Adj : so(3) so(3). 
For B £ sr)m(3), then using the cyclic property of the trace, we have 

(AAdsC)^- (AdsA C)^, 

i.e. Ads is self-adjoint with respect to the Frobenius inner product. Hence, 
9qT(0,/) is self-adjoint, and 

{^Qr{0,I).^i,p):^ = -\B2\{AdBf3,Adjfi)r 

It follows from this identity that if Ads and Adj have trivial kernel on so (3), then 
9qT(0,/) is invertible. Since J is diagonal, Jij — piSij, we have 

{Ad J p)ij = {pi - Pj)pij 

and hence 

\\Adjp,\\jr > min \p^ - pj \ \ \p\\jr 

It follows that if the Hessian Bij and the inertia tensor of B2 are both diago- 
nal, and such that Bij is invertible with distinct eigenvalues and Jij has distinct 
eigenvalues, then 9qT(0, /) is invertible. Recalling that the assumption that Bij 
and Jij are diagonal can be imposed without loss of generality, an application of 
the implicit function theorem proves the following proposition. 

Proposition 4.1. Assume that the Newtonian potential Vi of i{Bi) has a critical 
point O such that the Hessian of Vi at O is invertible and has distinct eigenvalues. 
Then for any given body B2 such that the inertia tensor J , given by |^.^[ ), has 
distinct eigenvalues, then for small A > 0, there is a homothetic motion Fx = 
X{Qx + p) such that Fi{B2) is equilibrated in the Newtonian potential of Bi, i.e. 
equation holds. 

Remark 4.1. The existence of reference configurations satisfying the assumptions 
of proposition pTT] was shown in [Sj section 5.2]. 
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Figure 1. Two-body configurations of tlie type constructed in 
TlieoremlO 

We can now apply tfie results of section [31 in particular corollary 13.31 to deduce 
the existence of a class of static, elastic two-body systems. Figure [1] illustrated the 
type of configurations which are covered by this result. 

Theorem 4.2. Assume that the Newtonian potential Vi ofi(Bi) has a critical point 
O, located in the complement ofi{Bi), such that the Hessian ofVi at O is invertible 
and has distinct eigenvalues, and let B2 be a given body such that the inertia tensor 
J of B2 has distinct eigenvalues. 

Then, for sufficiently small G, there is a static, self- gravitating elastic two-body 
system in Einstein gravity close to the configuration U F\{Bo2), with F a 

homothetic motion as constructed in proposition \4.1\ 

4.2. Axisymmetric bodies. In this section we consider configurations which have 
a rotational symmetry. We may without loss of generality assume that i(Bi) and 
1(^2) have have centers of mass at the origin O of the cartesian coordinate system 
on M^, and that the axis of symmetry of the bodies is aligned with the z axis. 
Thus letting (j) be defined by tani/i = y/x, we have that 9^ is a symmetry of the 
configuration, in the sense that dij,XBi =0,^ = 1,2. 

Suppose the Newtonian equilibrium condition (|2.ip holds. Let Z — Z(G, A) be 
the solution to the projected system, constructed using ProDOsition l2 . 1 1 and let N be 
the normalized force map as in section [321 For the present purpose it is convenient 
to take N as defined in terms of the material frame, i.e. 

As defined, N takes values in the dual of the space of Killing fields of (R|,(5). 
However, using the Euclidean geometry of R|, we may consider N as taking values 
in the space of Killing fields. We have [N, 9^] = 0. One easily checks that the only 
Killing fields which have vanishing Lie bracket with d(j, are linear combindations of 

and dz- Since is a symmetry of the body, it then follows that N is proportional 
to dz- We remark that this can be seen directly from the fact that due to the axi- 
symmetry of the body, any load must be along the z-axis. 

From the above discussion, it follows that the two components cannot be sep- 
arated in the z-direction. Thus, any axi-symmetric reference configuration with 
two components must have the property that one component is located "inside" 
the other, see figure [21 In this figure, each point corresponds to a circle, i.e. the 
bodies constructed are achieved by rotating the regions shown around the z-axis. 
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Figure 2 . Slices through configurations of the type considered in 
theorem 14.41 Each point in the figure corresponds to a circle. 

In each toroidal object is in equilibrium near, or in a toroidal cavity in a 

larger object. 

Thus, in order to solve the equation N{G, Z{G, A)) = 0, it suflices to consider 
the dz component of N. Further, the only motions we need to consider are those of 
the form z ^ z + a. 

The terms in the derivative 9yiN(0, /) which we need to consider are thus of the 
form 

If this quantity is non-zero, then in view of the remarks above, the argument which 
proves Proposition 13.11 proves 

Proposition 4.3. Assume that the reference body B C is in equilibrium in 
the sense that h2. 1 a\] holds. Let Z = Z{G,A) be the solution to VbiuB2^ = 
constructed in Provosition \2. 1[ and let N(G, A) be the normalized force map defined 
by fH). 

Suppose B is axisymmetric, with its axis of symmetry aligned with the z-axis, 
and with center of mass at the origin O. Further, assume that B is a disjoint union 
B = Bi U B2 with the property that 

is nonzero, where Vi is the potential ofi{Bi). Then there is an e > and a smooth 
map G 1-^ A, [0, e) ^ A with A(G) of the form z 1-^ z + a for a G M, such that 

N{G,A{G)) = 

Remark 4.2. Proposition \4. S\ gives an example of a situation where the normalized 
force map N{G,A) has degenerate derivative at (0,/), but where the symmetries of 
the situation still allow us to apply the same argument as in the non-degenerate 
situation. 
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Figure 3. The force generated on a body placed ofF-center in a 
spherical, elastic shell. 

Theorem 4.4. Suppose that the reference configuration satisfies the assumptions 
of proposition \^.3\ Then, for sufficiently small G there is a static, axisymmetric, 
self- gravitating two-body system in Einstein gravity, close to the configuration i{B). 

4.3. A body in a spherical shell. As mentioned in the introduction, the New- 
tonian potential inside a spherical shell consisting of a homogenous material is 
constant. Considering only rigid bodies in the Newtonian theory, it is thus possible 
to place a small body at rest at an arbitrary position inside the shell. Here we point 
out that if we consider instead elastic bodies, then this general construction is no 
longer possible. On the contrary, we argue here that the only static configuration 
of this type consists of a spherically symmetric body placed at the center of the 
shell. 

Consider an outer outer shell of radius R. Take coordinates such that the center 
of the ring is at r = 0. Place the small body at distance d from the center on the 
2;-axis and consider the system in Newtonian gravity. 

Let us consider linearized elasticity for this system. There linearized deformation 
of the inner body is zero because the force vanishes in the interior. Due to the 
principle of actio est reactio, the force of the inner body on the outer shell is 
equilibrated, and hence the linearized elasticity equation has a solution. 

The Newtonian potential generated by the small body is proportional to 

1 dz ^, 1 , 
- + ^3 

The leading order after the spherical field is a dipole field, i.e. a. i — 1 spherical 
harmonic. Figure [3] shows the upper half of the shell, with the small body. The 
deformation and the resulting force are indicated in the figure. 

The linearized deformation generated by V will therefore generate inside as dom- 
inating contribution a ^ = 1 spherical harmonic gravitational field 

5V = az 



STATIC MANY-BODY SYSTEMS 



15 



However, in this field the smaU body at d can never be equihbrated (recall that the 
1 = 1 part is the dominating contribution provided the radius of the outer shell is 
suffciently large.) Since the linearized problem has no solution, we cannot expect 
a solution to the full non-linear problem to exist. 
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